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ABSTRACT: We improve and generalize in several accounts the recent rigorous 
proof of convergence of delta expansion - order dependent mappings (variational per- 
turbation expansion) for the energy eigenvalues of anharmonic oscillator. For the 
single-well anharmonic oscillator the uniformity of convergence in g e [0, 00} is proven. 
The convergence proof is extended also to complex values of g lying on a wide domain 
of the Riemann surface of E(g). Via the scaling relation a la Symanzik, this proves 
the convergence of delta expansion for the double well in the strong coupling regime 
(where the standard perturbation series is non Borel summable), as well as for the 
complex "energy eigenvalues" in certain metastable potentials. Sufficient conditions 
for the convergence of delta expansion are summarized in the form of three theorems, 
which should apply to a wide class of quantum mechanical and higher dimensional 
field theoretic systems. 
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1 Introduction 



A recent rigorous proof of the convergence of scaled delta expansion (inspired 
by and developing further an earlier work by Duncan and Jones @), has provided 
for the first time a solid basis for understanding the success of the delta expansion - 
order dependent mappings (DE-ODM in the following). 

However the proof given there has several unsatisfactory features. First, the proof 
does not apply as it stands to the double well case {u 2 < 0) while numerical evidence 
suggests that (optimized) delta expansion converges for coupling constant above some 
critical value (i.e., the central barrier lower than a critical height). 

Also, delta expansion converges numerically well even in the g — ► oo limit || 
but the proof given in Ref. [Q] is valid for finite g only. (In other words, uniformity 
of the convergence in the whole g (> 0) axis was not proven). Even for finite g the 
absolute value of the remainder was numerically found to be smaller than our upper 
bound, suggesting that a tighter bound was possible to get, even if the bound found 
there was sufficient for proving the convergence for any finite g. 

The proof of |3J furthermore does not apply if the coupling constant g is taken 
to be complex. As the energy eigenvalues of the anharmonic oscillator and those of 
the double well are related by analytic continuation in g, it is interesting to know the 
domain in the complex g plane in which delta expansion converges. Seznec and Zinn- 
Justin || analyzed these problems and argued in a non rigorous fashion that delta 
expansion or order dependent mappings converges in some complex domain of g. 

These shortcomings of our previous proof have clearly to do with the fact that the 
detailed knowledge on the analytic property of the energy eigenvalue as a function of 
complex g, as explored by Bender and Wu ||, Loeffel and Martin ||, Simon 0, has 
not been fully exploited in the proof.[] 

In the present paper, we improve on these points and generalize our convergence 
proof in several accounts. First we generalize the allowed form of the reparametriza- 
tion (conformal transformation) of the coupling constant used to construct DE-ODM. 
Secondly, the DE-ODM for the anharmonic oscillator /double well system is proved 
to converge for complex coupling constant g in a domain on the Riemann surface of 
E(g), similar to the one suggested by Seznec and Zinn- Justin f3|. Also the proof is 

1 We thank J. Zinn- Justin for stressing this point to us. 
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considerably simplified with respect to our previous work. 

Furthermore, a general set of sufficient conditions for the convergence of DE-ODM 
will be given for a generic physical quantity considered as a function of the coupling 
constant and whose standard perturbative series is known. 

To set the general background for the present work, some known facts about the 
scaling relation a la Symanzik, the relation between the energy eigenvalues in the 
single well and those in the double well oscillator, and the general analytic structure 
of E(u 2 ,g) in g, will be briefly reviewed below. 

Consider the anharmonic oscillator with q 2M potential described by the Hamilto- 
nian, 

H ~^ + Y q + 2M q ■ (L1) 
By considering the scale transformation q — > £g, unitarily implemented on the Hilbert 
space of states (as first suggested by Symanzik), one finds ||, ||] a remarkable relation 
for the energy, considered as a function of uj and g, 

E{u\g) = C 2 E{^\e M+2 9) (1-2) 

valid in the analyticity region for any complex number £ 0). This relation plays a 
powerful role in the study of analytic properties of E(u 2 ,g). 

An example of Eq. flL2| ) is 

E K (u 2 ,g) = g y^E K (g-^V, 1) (1.3) 

valid for each level K. Combining this with the knowledge that the expansion in uj 2 
yields an analytic regular perturbation series, one proves H the (uniform) convergence 
of the large g expansion, 

oo 

E K (l,g)/g^ M +V = £ ^/(M+i)^ (L4) 

n=0 

for complex g such that 

\g\>9o>0 (1-5) 
(note that #o can always be increased so as to ensure the convergence at the border). 
Also, the choice £ = e m ^ 2 leads to a periodicity relation for energy, 

E(u J 2 ,e^ M+1 ^g) = -E(u J 2 ,g). (1.6) 
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Furthermore, by an appropriate choice of £ in Eq. (|1 . 2| ) it is possible to relate the 
energy eigenvalues in the single well potential to those in the double well potential. 
Indeed, the choice £ = e i7r//4 leads to 

= e-^ 2 E K {u 2 ,e^ M+1 >/ 2 g). (1.7) 

Analytic continuation of E(g) = E(l,g) and the structure of its singularities 
in complex g plane have been studied with different techniques and levels of rigor: 
Bender and Wu Q and Bender et al. [f?[] used in their pioneering works the WKB 
approximation, together with some general arguments; Loeffel and Martin || studied 
analytic continuation of the exact solutions of differential Schrodinger equation, and 
Simon || used the general theory of linear operators in Hilbert spaces. 

Some of their main results are: 

i) E(g) is analytic in the whole cut g plane || [|; 

ii) the Riemann surface has a global M+ 1-sheeted structure: a full rotation of angle 

2(M + 1)tt at fixed \g\ brings back to the original value of E{g) j|. [|; 

iii) E(g) has an infinite number of pairs of isolated square-root type branch points^ 

(Bender- Wu singularities) @] on the symmetric positions with respect to the 
phase ±z(M + l)7r/2 that accumulate towards g = (in some Riemann sheets 
other than the first), with asymptotic phase e ± i ( M + 1 ) 7r / 2 g g- a ^ these singu- 
larities in complex g plane the crossing of energy levels occurs. 

iv) A full closed tour of angle 2(M + in g plane starting from Ex(g) but this 

time by crossing the Bender- Wu cuts several times, takes one back to the energy 
eigenvalue of a different level, Ex+nig), where n depends on the number and 
the way the Bender- Wu cuts have been crossed [|J. 

v) For small enough g, E(g) is analytic in a sectorial domain 0: let < 9 < 

(M + l)7r/2 then there is a G such that E(g) is analytic in the subset 

[g\ < \g\ < G, \Avgg\ <6} (1.8) 

2 A rigorous proof of the detailed structure of higher Riemann sheets is still lacking (to our 
knowledge): Simon Q proved that the only isolated singularities of E(g) are algebraic branch points 
with no negative powers, but he could only assume the non existence of non-isolated singularities 
other than g = and of other pathologies. 
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of the Riemann surface. Also in any such sector the Rayleigh-Schrodinger per- 
turbation theory is asymptotic uniformly in the angle H. 



vi) A unique strong coupling expansion Eq. (|1.4j) exists which converges uniformly 
for \g\ large ||. 

The rest of the paper is organized as follows. We first recall in Sec. |2] the relation 
between the delta expansion and the order dependent mappings, which serves also to 
review briefly these ideas and at the same time to define our convention. The con- 
vergence proofs are given in Sec. ^, which constitutes the main result of this paper. 
In Sec. H convergence of the optimized expansion is discussed. The case of general 
q 2M oscillators is discussed in Sec. |[ Sec. || discusses briefly the convergence of delta 
expansion for the Green's functions of the anharmonic oscillator. Summary and Dis- 
cussion is given in Sec. [7[ Several technical issues and examples are collected in Ap- 
pendices: Appendix A| discusses some properties of the conformal transformation used 
in the proof, |Appendix B| the convergence proof of DE-ODM in a zero-dimensional 



analogue model. The condition on the trial parameter following from the optimiza- 
tion condition is discussed in [Appendix Q . In Appendix D| we present the calculation 



of critical exponents in 0| theory by use of DE-ODM. Finally a reconstruction of 
the strong coupling expansion coefficients from the standard perturbation series is 
discussed in Appendix E . 



2 Delta expansion versus order dependent mapping 

The delta expansion is a method to perform a systematic calculation, taking as 
the zeroth approximation an unperturbed action v4 ({fi}) which contains a set of trial 
"vacuum" parameters {Q} (e.g., trial frequency, mass)|] Perturbation with respect 
to the rest of the action 5[A — Aq({Q})] is calculated in the standard fashion, and 
the iVth order result (with 5 = 1) S^ E ^ is computed by fixing {Q} order by order by 
some appropriate criterion, e.g., by an optimization procedure. 

For the quartic quantum mechanical oscillator Eq. (|1 . 1|) with M = 2, the "unper- 

3 In this respect the delta expansion is close in spirit to the Hartree or self consistent approxi- 
mation. There are many works related to delta expansion. Standard earlier references can be found 
in 0]. See B for more recent ones. 
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turbed" and "interaction" Hamiltonians can be taken as 

r? tt 2 uj 2 -Q? a 

H (n) = ^ + '-L q 2 ; Hj(n) = ^-q 2 + ^ (2.1) 

In this case for the energy eigenvalues was proved [[!]] to converge to the exact 

answer as iV — > oo if the trial frequency is scaled with the order iV as 

XN = ^L = CN\ (2.2) 

UJ 



where either 



1/3 < 7 < 1/2, C > 0, (2.3) 



or 



7 = 1/3, C > a c g 1/3 , a c = 0.5708751028937741 • • • . (2.4) 

Also the divergence of delta expansion with xn lying outside the above range, was 
explicitly shown. The convergence of the "optimized" delta expansion (in which 
VL N is determined for instance by the condition dS^/dO.^ = 0, at each order), is a 
consequence of the above (see Sec. f| below). 

Formally, the delta expansion in this case can be generated by a substitution, 

uj - n ( 1 + 5 ^f-j = u,x[l - 5(3(x)] 1/2 , 
9 9, 

x = -, P(x) = l-\, (2.5) 
UJ x z 

in the standard perturbation series 



OJ 3 



E^\u\g) = ujY^Cu 

n=0 

= ujE^\1, 9 /uj% (2.6) 

followed by the expansion in 5 (and 8 — 1; u — 1 at the end). 

In the order dependent mapping method of Seznec and Zinn- Justin ||, one con- 
siders a physical quantity S(g) regarded as the function of a (in general complex) 
coupling constant g, and whose perturbative series in powers of g is known up to the 
A^th order. One makes then a conformal transformation, 

9 = pF(0) (2.7) 
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to the new (in general complex) variable (3. The (order by order-) adjustable param- 
eter of the transformation, p, is kept always real and positive. The transformation is 
such that g = and g = oo are mapped to (3 = and (3 = 1 respectively; the full 
range [0, oo] in g is mapped to [0, 1] in the space of (3. For instance, in the case of the 
anharmonic oscillator the choice 

m - ir J ^. (2.8) 

with a > a fixed parameter, will work. The iVth order approximation S^ DM ^ is 
obtained by re-expanding S(g) in powers of (3 up to iVth order, 

N 

S(g)^S { ° BM) = £P„(p)/r, (2.9) 

n=0 

and p is fixed by some criterion, order by order. In Ref. || Seznec and Zinn- Justin 
adopt the criterion of the fastest apparent convergence, i.e., they require p at order 
N to be a root of Pjv+i(p) = 0. 

In the case of the energy eigenvalues of the quartic anharmonic oscillator, the 
same authors considered the function ty(/3,p) related to the energy by 



1 
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E ^9) = Q-pyfi [^-P) 1/2 E(l,g((3))\ = *(p,p); (2.10) 

with a particular conformal transformation 

^ ) = /? (T^)^ ; p = ^wr (2 - n) 

Only the function \l/(/3, p) (but not the prefactor 1/(1 — (3) 1 / 2 = x) is re-expanded in 
(3. 

Comparison between the set of equations ( |2.5|) and ( |2.6| ), and the set Eqs. ( |2.10| ) 
and ( P-ll|) , shows that in this case the delta expansion and order dependent mapping 
(with the particular conformal transformation chosen above) coincide formally. (3 in 
the two cases can be identified. 

It is clear that the two approaches are closely related to each other. The original 
form of the delta expansion has a certain intuitive appeal, being physically motivated. 
However it is not obvious how to apply such a method to higher dimensional field 
theory models where mass renormalization is important. The order dependent map- 
pings, on the other hand, being mainly mathematically motivated, appear to be more 
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flexible. In particular the conformal transformation in the latter can be of a more 
general form than suggested by the shift of the quadratic term. In the present paper 
these methods will be referred to indistinguishably as DE-ODM. 



3 Convergence proof 

In this section convergence proofs of DE-ODM are given. The convergence in the 
case of the quartic potential (M = 2) for finite positive g is given in Sec. |3.1| , for a 
wide class of conformal transformations, Eqs. ( ^.lp - flOp . In Sec. |3.2| , the convergence 
proof for quartic potential is generalized to complex values of g (including higher 
Riemann sheets of E(g)) by use of a particular conformal transformation (a = 3/2 in 
Eq. (|3.1| )). All the (sufficient) conditions used for the convergence proof of DE-ODM 
are summarized in Sec. [O] in the form of three theorems to be applied for a generic 
function S(g) of g. 

The strategy for the proof of convergence is common in all cases, and can be 
summarized as follows (with reference to the energy of anharmonic oscillator). 

1) Consider the energy as a function of the new variable j3, E(g(/3)), through the 
conformal transformation, 

9 = pF(f3) = p ^ \ (p real positive). (3.1) 

For consistency with delta expansion notation (see Sec. it w iU be useful to 
define also the variable x by: 

The parameter a will be always chosen in the range 

1 < a < 2 (3.3) 



to ensure some properties essential for the proof (see |Appendix A|) . Since E(g) 



is analytic in the (open) cut g plane (see Sec. [I]), the composed function E(g(f3)) 
is analytic (at least) in its (open) image D in j3 plane. It is shown in |Appendix A| 



that T> is the interior of an apple like figure minus the negative real axis. (See 
Fig. 1 for the case a = 3/2.) 
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2) Introduce a function \I/(/3, p) by 

E W)) = 7T30V*GM- ( 3 - 4 ) 

-£?(<?) in the case of the energy of g 2M -oscillator is known to behave as ^ 1 /( M+1 ) 
at g — > oo (see Sec. [TJ); the extraction of the factor 1/(1 — is made so that 
the function ty(/3,p) is bounded at /3 = 1. The condition for that is: 

* > TT^-T- (3-5) 

In cases we are interested in the convergence proof for values of g in higher 
Riemann sheets, more special choices of a will be required; see below. 

By construction p) is analytic in T>. 

3) Write Cauchy's theorem in the (3 plane (See Fig. 1), 

where T is a closed contour which traces the boundary of the domain T> (we 
assume that the energy is continuous on the two edges of the g cut, see below). 
Note also that the contribution to the Cauchy integral from the portion of T 
that turns around (3 = is negligible because the uniform asymptoticity of the 
perturbative series for E(g) implies that E(g) —>■ 1/2 for \g\ — ► 0, see Sec. [1]). 
In certain cases one can (and must) enlarge further the contour in the (3 plane 
(depending on the analytic structure of p)); such a deformation is crucial 
in the proof of convergence in higher sheets. 

4) Taylor expand ^(f3,p) in (3 up to the iVth ordeifj 

N j 



n=0 n[ 




(3.7) 



4 Derivatives respect to (3 in (3 = are taken from the positive axis and can be exchanged with 
the integral in /?' of Cauchy formula essentially due to regularity of ^ in this region (in particular 



the behaviour of Disc "J for /3' — > , see Sec. 3.1) 
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The remainder for the energy, Rn{/3) = [^(P,p) — ^n{P, p)]/(l — PY is then 
given by 

p N+l r dp' *(/3',p) 



Rn{(3) ~ & 2™ - P) ■ m 

Study of various contributions (essentially from the two regions /?' ~ 1 and 
P' ~ 0) in Eq. ( |5TS| ) in the limit iV — >• oo, p = g/(x 2a P) — > 0, leads to convergence 
proofs and the related conditions for various parameters. 



3.1 Proof for the quartic oscillator (M = 2) for real g 

We first prove the convergence of DE-ODM to the energy eigenvalues of the quartic 
oscillator (M = 2 in Eq. ( |1 . 1| ) ) for any finite positive value of g. Following the general 
steps outlined above, we make the conformal transformation 

g = p F(p) = p _f^ a i (P real positive). (3.9) 

with 1 < a < 2. Finite positive real g will be mapped in P <G (0, 1) by this conformal 
transformation. 

We introduce the function ^{P,p) by ^(P, p) = (1 — P) a E(g(P)) (with a restricted 
as in Eq. ( |3.5| )) and use the Cauchy's theorem with the contour given in Fig. 1. 

To analyze the remainder 

we divide the contour on f3 plane in two parts: Ta the part which wraps around the 
subset [P c = — l/(a — 1), 0] of the negative real axis, and T B = T/Ta the rest of the 
contour, the upper and lower arms of the apple like curves in Fig. 1 (a complete study 
of T for generic a is given in |Appendix A| ) . 



Consider first the contribution from Fa- Clearly the factor (1 — P') a in 

p) = (1-PTE (g(P')) (3.11) 

takes the same value on the upper and lower sides of the cut along j3 G [P c , 0]. Hence 
the discontinuity (imaginary part) of ^ can be expressed in terms of that of E and 
this part of the remainder can be written as 

R (A) p n+i rdp> (i-pri™E(g(p')) 
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In the limit N — > oo, p — ► 0, we have 

g(P') = p- 



1-/3') 



for any value of f3' G Ta- Thus the known g — > behavior of lm.E(g) 



~ C E (-g) c exp 



(3.13) 



(3.14) 



(-^) 1/b . 

can be used in the integrand (b = M — l = l,a = 4/3 here and c is a positive constant 
depending on the energy level considered). One then finds 



\Rn I — cons t- 



IJV 



/oZ-CJV+l+c) 



exp 



a 

w 



(3.15) 



|l-0|* r 

where the constant (const.) is independent of the order of DE-ODM and of g,f\ and 
inequalities 

|0 - /3'| > |0| (3.16) 



and 



|l-0d > |1-01 > 1, 



1-Pc 



< W)\<pW\ 



(3.17) 



valid for 0' G T^, G [0, 1] have been used. Taking an upper bound by raising the 
upper limit of integration to infinity one gets: 



I l 

\ K N I 



IJV 



< const. 



P 



N+c 



T(N + c) 



(pN) 



N 



N 



X 



2a 



N 



where use was made of the asymptotic relation, 



0=1-4-1- r- 



l/a 



(3.18) 
(3.19) 

(3.20) 



Next consider the contribution from Tb- The range of a, 1 < a < 2, has been 
chosen so that the following inequalities 



|0 -01 > 



(3.21) 
(3.22) 



Any such constant will be indicated this way in the following. 
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hold for f3' on V b (see |Appendix A| ) and (3 G [0, 1]. 

Furthermore we assume the continuity of E(g) on the two edges of the negative 
g cut. This physically reasonable assumption is motivated by the fact that all the 
known singularities of E(g) in the nearby domain 



g = \g\ e 



7T < 9 < 



3tt 



or < 9 < -7r 

2 ~ ~ 



(3.23) 



are simple square- root type branch points [§]] (see footnote 0). Also, due to the fact 
that the energy is analytic on the cut plane, that E(g) ~ g 1 ^ 3 for large g, and that 
E(g) is bounded in the limit g — > 0, the following bound holds on the cut plane 
including the two edges of the negative <? axis: 



1^)1 < (X + Slpl 1 / 3 ). 
In terms of \I> we have (on the whole of F U £>): 

Pl/3| 



(3.24) 



l*GMI< 



1/3 



A + B 



|l-/?l c 



ii-/?r 



(3.25) 



from which (by use of Eq. (3J5)) one easily derives 

< const.; (3'eT B . 



(3.26) 



It follows from Eqs. fl3.21| ), ( |3.22| ) and ( |3.26[ ) that the contribution from to 
Rjsr, R N , is bounded by 



It\ U) \ < const. J^T?. ~ e"^ 2 ~ exp 



l+o- 



-AT 



l/a 



(3.27) 



(where Eq. (|3.20|) has been used). 

It clear from Eqs. ( gl9|) and (|3~27|) that both pieces of the remainder, and 
R N , vanish in the limit of large N if x or p scale with N as 



x = CW- ^<1<\ (C>0) 



or 



iV7' 



1< i < a (C > 0). 



(3.28) 
(3.29) 
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Accordingly DE-ODM for an energy eigenvalue for the quartic oscillator, constructed 
with the class of conformal transformations in Eq. ( |3.9|) with 

1 < a < 2; (3.30) 

and 

<t>~; (3.31) 

converges pointwise to the exact answer for any positive g, if x (or p) is scaled as 
above. 

Finally, for completeness we state here that the convergence also holds if x (or p) 
is scaled precisely as 

x = CN i/(2a). or p = ^-; (3.32) 

with 

C >(^) or C'<±. (3.33) 



This is a special case of the Theorem 3 treated in Sec. |3.3| and is a solution 



of the set of equations (p.77|) and (|3.80| ) (with b = 1). In this case convergence is 



ensured by the factor \/3\ N appearing in both pieces of the remainder. For a particular 
conformal transformation with a = 3/2, it reduces to the condition with (Su^/A) 1 ^ 3 = 
0.5708751089 • • •, found earlier 0. 

Remarks 

(i) The result of this subsection generalizes considerably the result found earlier (the 
case with a = 3/2) by the present authors (]J and argued for by others ||, 
allowing a class of conformal transformations to be used in the construction of 
DE-ODM. The convergence proof is also simpler than our original one. 



(ii) As is suggested from the bounds Eq. (|3.19[) and Eq. (|3.27|) , the convergence proof 



of this section can be actually extended to any complex g in the first Riemann 
sheet satisfying Rel/g 1//a > 0. See Theorem 1 in Sec. pT3| . 



3.2 Proof for the anharmonic oscillator/double well system 
(M = 2) with complex g 

We consider the quartic oscillator (M = 2) again but now the coupling constant 
g is allowed to lie anywhere in the Riemann surface of E(g). Through the scaling 



12 



relation Eqs. (|1.2j) and (|1.7| ) this permits us to study simultaneously both the single 
and double well problems. 

As before we make the conformal transformation, 

g = p F(f3) = p _f^ a i (P real positive). (3.34) 
and introduce the function \I/ by 

E W)) = k=0F*(P>p)- ^ 3 - 35 ) 

The key observation which permits us to extend the proof to g lying on higher 
Riemann sheets, is the following. The existence of the large g expansion for E(g), 

oo 

E K {l,g) = g L/ *'Ed n (g- 2 '*) n , (3.36) 

71=0 

with a convergence radius go, shows that for the following particular choice of param- 
eters 

a = |; a = -k, (A; = 1,2,3,---) (3.37) 

the function \l/(/3, p) is analytic^] in (3 inside a certain little circle C of radius r centered 
at /3 = 1, as well as in the domain T>. The first of the conditions in Eq. ( |3.37| ) undoes 
the branch point at g = oo due to the overall gf 1 /( Af + 1 ) factor; the second eliminates the 
same branch point due to the fractional powers of g in the expansion parameter |3J. 

Note that such elimination of the branch point at /3 = 1 is possible only for a par- 
ticular choice of the conformal transformation. Indeed, the condition for analyticity 
in C, Eq. ( |3.37| ), together with the condition on the exponent a of the transformation, 
1 < a < 2 (which is needed to keep (3' along Tb outside the unit circle), implies that 
the only possible value for k in Eq. ( |3.37|) is A; = 1. The only allowed conformal 
transformation is the one with 

a = -. (3.38) 

From the first of Eq. ( |3.37| ) it also follows that a = 1/2. These particular values are 
assumed hereafter in this subsection. 

The radius r$ of C is easily found to be, in the limit iV — > oo, 

«~ £ , (3.39) 



By use of Weierstrass' theorem on series of uniformly convergent analytic functions. 
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by use of Eq. Q). 

Since now the function ^((3, p) is analytic in the joint region T> U C, the contour T 
can be enlarged to T' so as to trace its boundary. See Fig. 2. As will be seen shortly, 
this apparently minor deformation of the integration contour is crucial. A further 
enlargement of the radius of C is not possible because the images of the Bender- Wu 
branch points appear just outside C.Q The smallest distance from one of those to the 
point (3 = 1 shrinks to zero together with r$ as iV — > oo. 

Before proceeding further, some subtleties involved in analytic structure of the 
function \1/ and in the consequent enlargement of the contour r — > V should be ex- 
plained. First let us note that conformal transformation g — p(3/(l — P) 3 ^ 2 introduces 
the doubling of the number of the global Riemann sheets (six) in E(g(j3)) considered 
as a function of (3 as compared to the situation for E(g) (three). As is clear from 
the form of the strong coupling expansion Eq. ( |3.36|) , there are a square- root branch 
point at (3 = 1 and a (global) cubic-root branch point at (3 = {(3 = oo is a sixth 
root type branch point). 

A detailed study of the conformal mapping from the g Riemann surface to (3 
Riemann surface shows the following correspondence. Two values of (3 with the same 
(geometrical) position on a pair of (3 Riemann sheets, connected by going around 
(3 = 1 once (which may be called "mirror" Riemann sheets), are mapped to values of 
g differing by the phase e 3?r \ For example, the domain T> of Fig. 2 corresponding to the 
first and its mirror (second) sheet, is the image of the cut plane — tt < Argg < tt and 
the angular region 2tt < Arg g < Att, respectively. Analogously, the part of the circle 
C which lies outside T> (Fig. 2), represents the image of the region tt < Axgg < 2tt; 
\g\ > go and that of 4n < Argg < 5n = —n; \g\ > go- 
Now, the analyticity of ^ in C around (3 = 1 implies, due to the uniqueness of ana- 
lytic continuation, that the function \l/ takes exactly the same value at corresponding 
points in each pair of mirror Riemann sheets. 

In terms of E = (1 — (3) l / 2 ^> this implies a relation 

E(e 3m g) = -E(g) (3.40) 

for any g. (This relation of course follows immediately from the scaling relation, 

7 <7o is larger than the largest absolute value of the Bender- Wu branch points. For otherwise there 
would be values of g such that \g\ > go and for which the strong coupling expansion for E/g 1 / 3 does 
not converge, which is a contradiction. 
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see Eq. (|1.6| ). Nonetheless the above derivation, which essentially hinges upon the 
existence of a strong coupling expansion, can be used in a more general situation: See 
Sec. p73| .) 

The true meaning of "enlarging the contour from T to T'" can be appreciated now: 
As j3' moves along V twice, first in the first sheet and continuing in the second sheet 
and then back, its image g' goes through the full g Riemann surface once, as shown 
in Fig. 3. 

Nevertheless, in applying Cauchy's theorem for ^((3, p) one may ignore the double 
looping of the V contour: on whichever sheet a given f3 may lie, the integration over 
one of the loops vanishes, the factor — f3) being analytic inside it. 

Coming back to the convergence proof, for any (3 lying inside T> U C we study the 
remainder 

P N+1 r dp 



Rn{(3) ~ (l- & 2ri{pr +i (P'-py (3 ' 41) 

By splitting again the contour in two parts, Ta the part wrapping around the negative 
real axis (common to T and T'), and T' B = T'/Ta, the rest of the contour (see Fig. 2). 



The contribution from Ta can be treated as done in Sec. [O, except that a = 3/2 



and variables g and (3 are complex now. By use of Eq. ( |3.17| ) and of 



\(3-f3'\ > Kp = max(|Im/3|,Re/3), (3.42) 

one finds easily 

\B\ N+1 / n\ N+c 

\R { ^\ < const. - p-c(tL) T{(N + c)) 



Kp\l-f3\y 2P \aj 
- ^^~\9\ 1/3 (PN) N , (3.43) 

where an asymptotic relation (3 ~ 1 — (p/g) 2 ^ 3 has been used. 

The appearance of K@ in the denominator of the first line of Eq. ( |3.43| ) means that 
the values of (3 in T> U C which are arbitrarily close to the real negative axis, must 
be excluded. We shall see shortly that the remainder vanishes as iV — > oo only for g 
lying on a certain region of the Riemann surface (Eq. ( p.49| )). For any g there and for 



any N > N (with iVo independent of g) the inequality Kp > const. \(3\ holds, hence 
the convergence will indeed be uniform in g in that region. 
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( B) 

There is an important difference in the estimate of R N due to the enlargement 
of the contour from r to V. For f3' on T' B , the following inequality holds: 

|/3'| >l + r o cos(0 a ) + O(r o 2 ) (3.44) 

where cb a = tt(1 — 1/a) is the angle between the apple like curve and the positive 
(3 axis in (3 = 1, see |Appendix A| . For the generic range 1 < a < 2 of a, the angle 
(f) a lies in the region < <p a < it/ 2 so that cos^q, > 0. For the particular conformal 
transformation with a = 3/2 we are considering here, (b a = 7r/3; hence cos<b a = 1/2. 

Also, to bound the factor l/\/3' — /3|, one must restrict g to range in a complex 
domain such that the corresponding variable (3 satisfies 

|/3'-/3|>er V/3' eT' B , (3.45) 

where e is a small positive fixed ((3 independent) constant of order 0(1). Appearance 
of Tq (radius of C) in the condition above is due to the fact that for any fixed finite g 
its image (3 approaches 1 in the limit N — > oo and that the nearest f3' are those on the 
boundary of C. In terms of g, it is sufficient that g is restricted to be separated from 
the inverse image of the contour T' B (Fig. 3), by a fixed (g independent) finite distance. 
It will be seen below that the condition for the (uniform) exponential vanishing of 
the remainder automatically ensures such a condition to be satisfied. 

Furthermore, since E(g) ~ g 1 ^ 3 at large g on the whole Riemann surface, ^ is 
bounded in the neighborhood of (3 = 1 (see Eq. (|3.5| )): 



W,p)| < const.; (3'eT B . (3.46) 
Collecting these results one gets 



JV+l / 1 \ N + 1 



\R N I < const. 



11/3 




exp -iVp^ Re— - + — - , (3.47) 



where Eq. (|3.39| ) has been used. Note the presence of the (^-independent exponential 
factor in Eq. (|3.47|) which was absent in Eq. (|3.27| ). 

From Eqs. ( |3.43| ) and ( |3.47| ) one sees that the full remainder R^' + R^ vanishes 
(uniformly) in the limit iV — > oo with 

p = -^ with l<y<3/2, (3.48) 
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if the complex values of g are restricted in (any closed subset of) the region 

Re + — U > °- ( 3 - 49 ) 

9 2/3 2g 2/S 

The regions excluded by such a condition are one heart-like closed domain 7i\ 

io/q 2/^ 2$ 3tt „ 9tt 

\g\ 2/3 <-2g 2 /3 cos—; — < 9 < —■ 9 = Axgg, (3.50) 

and an analogous domain 7i 2 with 9 — > —9. See Fig. 3. The boundary curves of 
these domains cross the straight lines corresponding to Argg = tt, 2tt, An (= — 2tt), 
5n (= — 7r) at modulus \g\ = go. 

As is clear from Fig. 3, the (inverse) image of the contour T' lies inside Tii U H.2- 
It follows that the condition for the uniform vanishing of the remainder, Eq. (|3.49 ), 
automatically guarantees Eq. ( 13.451) . 

We have thus proved that DE-ODM for an energy eigenvalue of the quartic os- 
cillator, constructed with the conformal mapping with a = 3/2 and with o = 1/2, 
uniformly converges to the correct analytic continuation E(g) in every closed subset 
of the region 

Re 4n + — U >0 > ( 3 - 51 ) 

if p or x is scaled as 

p = ^; K V < |, (3-52) 

or 

x = cg l/3 N^, - < 7 < - (3.53) 
3 2 

(C or c is real and independent of g and N). 

As in the previous case, the convergence holds for the extreme case of scaling with 



a 

p = — ; or x = eg 



i/3jyi/3 



c>c* = 0.5708751089 • ■ ■ or C < 4- (3.54) 

Remarks 

(i) The g independent factor in the exponent of Eq. ( |3.47D shows that even for real 
finite g, the bound for the remainder is actually numerically smaller than sug- 
gested by Eq. (|3.27 ), although the exponential behavior (\R\ < exp(— const. iV 1-27 ' ' 1 
is the same. The convergence is fastest for the choice 7' = 1. 
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(ii) The same g independent factor is especially important for the strong coupling 
limit, g — > oo. By using Eqs. (|3.43|) and (|3.47|) we prove that DE-ODM con- 
verges to the correct answer for lim^oo E(g)/g 1 ^ 3 , if x or p is scaled according 
to Eq. ( |3.53| ) or to Eq. ( pj.54 ), and if the limit g — > oo is taken at each order in 



N, before N is sent to oo. This remarkable fact was noted earlier empirically. 

(iii) As is noted in Introduction, E(g) continued analytically to g = \g\e 3m / 2 is (apart 
from a phase factor —i) equal to the energy eigenvalue of a double well. The 
result above proves^ that DE-ODM for the quantum double well converges to the 
correct answer as long as \g\ > 2 3//2 go- This explains the empirical fact that lower 
order (optimized) delta expansion for the double well gives good approximate 
results for coupling constant above some critical value (or the central barrier 
below a critical height) but fails below it .0 From our own numerical study 
of DE-ODM for the double well up to iV ~ 50, we find a numerical value 
go — 0.10 0.15. Due to the fact that Eq. ( 3.51|) refers to an upper bound 



of the remainder, this should be considered as an order of magnitude estimate 
for the convergence radius of the strong coupling expansion. See the following 
remark. 

iv) Our result depends on the knowledge of analyticity of E(g) in the cut g plane: 
Bender- Wu singularities are confined in the angular region, 7r < Argg < 2tt; 
\g\ < go (and a similar region with g — > —g), but no further details on their 
location have been used. A tighter bound (and a wider convergence region 
in g) can be obtained by a further deformation of the contour near (3 = 1, 
if E(g) turns out to be actually analytic in a wider region. Indeed, if E(g) 
would be analytic in the region, |Argg| <7r + A, 0<A< tt/2, then l/(2g^ 3 ) 
in Eqs. (|3~4p , (|3~49|) , (|3T5T|) would be replaced by cos[(tt - 2A)/3]/^ /3 - The 
estimate of go from our numerical study of delta expansion for the double well, 
would accordingly be changed to 0.10 < g < 0.42, which seems to be consistent 
with the known strong coupling expansion coefficients ][[|, 11]. 



8 The Nth order delta expansion approximant for the double well is precisely |jj given by the 
substitution x 2 — > ~x 2 (that is, oj 2 — ► — u 2 ) made in the corresponding approximant of the single- 
well anharmonic oscillator. It follows that the remainder for the double well is correctly given by 
Rjf + Rft apart from an overall factor i. 

9 This was noted independently by Okopinska by Kleinert |fjo|| , by the present authors (un- 
published), and by E. Hofstetter (private communication). 
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(v) E(e lw g) with real positive g corresponds to the energy in an unstable potential 
well. According to our result, the complex energy "eigenvalues" of such a system 
can be correctly reconstructed by delta expansion/order dependent mapping, as 
long as g > go (sliding regime, rather than tunnelling). This nicely explains the 



empirical success of a variational-perturbative approach discussed in Ref. fll3 



3.3 General sufficient conditions for convergence 

In this subsection the (sufficient) conditions for the convergence of DE-ODM are 
summarized in the form of three theorems. 

Theorem 1: Let a function S(g) be given such that: 

1) S(g) is analytic in the complex g plane cut along negative axis, continuous on the 

two edges of the cut and bounded at g = (on this sheet); 

2) S(g) ~ g p with p > for g — > oo on the cut plane; 

3) The discontinuity of S(g) along the cut behaves as 



DiscS(g)~Cs(-g)- c exp 
for g — > 0~ with a > 0, 1 < b < 2, c real, 



(-g) 1/b 



(3.55) 



then the sequence {S^lg)} of approximants for S(g), constructed with conformal 
transformation 

g = p F(0) = p — f ^ real positive) (3.56) 

with 

1 < a < 2, (3.57) 
followed by a Taylor expansion in (3 up to order N of 

*{P,p) = (\-p)*SW)) (3.58) 

with 

a > ap, (3.59) 
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converges to S(g) as N —>■ oo uniformly in every compact subset of the region 

Re^>0, (3.60) 

gl/a V I 

if the real positive parameter p is chosen to scale with N as 

p = -^ 7 ; with b<i<a. (3.61) 

(C independent of g.) 

Theorem 2: Let S(g) be given which satisfies l)-3) of Theorem 1 and moreover: 

2') its large g expansion 

oo 

S(g) = g P Y, d n(9- q ) n (3.62) 

n=0 

converges uniformly for \g\ > go and there exists a strictly positive integer k 
such that 

q < k <2q. (3.63) 

Then the sequence {5jv(g)/g p } with Stf(g) constructed as in Theorem 1 with the 
conformal transformation Eq. (\3. 5dJ with a = k/q (thus lying in the range 1 < a < 
2), by use of the function ty((3,p) with o = ap, converges to S(g)/g p as N — > oo, 
uniformly in each closed subset of the domain (including g = oo) 

Re + 4^ cos - i / a ^ > ° ( 3 - 64 ) 

9 g ' 



for any choice of scaling Eq. ( \3. 61\ ) of the positive parameter p = C'/N 7 ' . (C inde- 
pendent of g.) 



Proof of Theorem 1: To prove Theorem 1 one proceeds as in Sec. |3" H3.1| : using 



1) and 2) one can write a dispersion relation as in Eq. ( |3.6j ) for ^ from which the 
expression Eq. (|3.8f) for the remainder follows. In particular from 1) and 2) one 
derives \S(g)\ < A + B\g\ p in the cut plane including the edges of the cut along the 
negative real axis, hence (with the choice Eq. ( |3.59| )) ^ is bounded by a constant on 
r U T> (r and T> are defined in Sec. |3|). The integration contour T is then split in 



and as in Sec. 3.1 
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As for Rpp, one expresses the discontinuity of \I/ along (3' G [/3 C , 0] by that of 
S(g(/3')), obtaining 

M) (3 N+1 fo d(3> (1 - /?') CT Disc S(g(P')) 



Rn ~ (1 - 4 27rz - 0) ' (3 ' 65) 

Then by using Eq. ( |3.55| ) for Disc S(g) (remember that g{(3') — > 0~ as A" — > oo for 
(3' G [ftc, 0]), the inequality 

\P-P'\ > const. |/3|, (3.66) 
(which is valid for g in ( |3.60| )) and the fact that the inequalities 

1 < |1 -fi\ < const., 

const. p\p\ < \g{(3')\ < p\(3'\ (3.67) 

hold (uniformly) on Ta, one finds an upper bound as was done in Eq. ( |3.15| ). Ex- 
tending the integration range up to /3' — oo one gets (on the chosen compact): 

, ks \P>\ N / n \ N+c 

[iffll < const. , p~ c (4) r(6(AT + c)) ~ M^piVT- ( 3 - 68 ) 



|1-/?K Va fe , 

As for contribution we observe that the choice of a Eq. ( p.57| ) is essential to 



derive the inequality (see |Appendix A| ) 



\P'\ > 1. (3.69) 
Also, for each compact subset of the domain Eq. (|3.60| ) the inequality 



\/3 - (3'\ > const. \l-/3\ (3.70) 
holds (uniformly in (3, (3'). One immediately gets (exploiting the boundedness of \I/) 

\R\i N + l ) 



\Rn ] \< const. jY^gjl^ ~ ex P 



-NRe - 
\9 



(3.71) 



Theorem 1 easily follows from Eqs. (|3.68| ), (|3.71| ) and ( |3.61| ) (note that the upper 
limit on 7' comes from Flffl while the lower limit from R^). Q.E.D. 

Proof of Theorem 2 The new hypothesis 2') on the existence of a large g 
expansion of S(g), together with the choices aq = k and a = ap, ensures (see 
footnote 0) the analyticity of \I> in the extended region DUC enclosed by T'. 
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As in Sec. [3.2| , the analyticity on the finite disk C around (3 = 1 together with the 
uniqueness of analytic continuation implies that the function ty((3,p) takes the same 
values on all the "mirror" Riemann sheets, i.e., those connected by going around the 
branch point at (3 = 1. As a result, on applying Cauchy's theorem one can work on a 
single sheet corresponding to the given value of g; loops on all other Riemann sheets 
give vanishing contribution, enclosing no poles. 

We proceed as in Sec. [3.2| , by splitting T' in and T' B . The treatment of 
contribution is the same as in Theorem 1 and one finds Eq. ( |3.68|) . 

As for the contribution of T' B , we observe that the condition Eq. ( |3.63|) (enforcing 
1 < a < 2) and the chosen range of g, Eq. ( p.64j ) , lead to the inequalities 

\P' ~(3\> er Q - 

|/3l>l + r o cos(0 Q )+O(r o 2 ) 
for (3' G T' B . (e is a positive (3, (3' independent fixed constant.) Together with the 



boundedness of \1/ (which follows from 2')), these inequalities give, as in Sec. \5?2 



\Rn I — cons t- 



1 — (3\ a r \ 1 + r cos i 



M P exp 



-Np 1/a ( Re -^r- + cos < 

V 9 l/a gl ,a 



(3.72) 



The theorem follows from Eq. (|3~68|) and Eq. (|3~72|) . Q.E.D. 



The convergence proof expressed in the two precedent theorems actually holds 
if the trial parameter p (or x) is scaled with the order N precisely as p ~ (or 
x ~ A^ 1 ^ 2 ")), under certain conditions on the proportionality coefficient. Since this 
case was discussed earlier by several authors in connection with the optimized DE- 
ODM for the anharmonic oscillator, we state the result for this extreme case as a 
separate theorem.^ 

Theorem 3: Theorem 1 and Theorem 2 hold if (maintaining the other respective 



10 Note however that in general such a particular scaling of p or x does not follow from the opti- 
mization condition alone (either principle of minimum sensitivity or fastest apparent convergence) 
used frequently in the literature. It is only after imposing a further condition, e.g., of choosing the 
root of smallest modulus that optimization equation fixes the particular scaling of x. See the related 
discussion in Sec. |]. 
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hypotheses) we scale the parameter p, instead of Eq. (\3. 61 ) as 

a 



< C < ( — ) (3.74) 



and 

< C < ( 

where u* is a function of a, b defined in Eq. ( \3. 8(\ ) below (see also Eq. ( \3. 82J for less 
stricter but a, b independent range of convergence) and a, b are those of Eq. ( \3. 53j ). 

Proof: We proceed exactly as in Theorem 1 and Theorem 2, splitting the re- 
mainder in two pieces and R^ . With the scaling Eq. ( $773$ ) the bounds on R^ ] 
Eqs. ( |3.71| ) and ( |3.72|) still hold unmodified and enforce convergence of that contribu- 
tion in the corresponding range of g. It will be shown that with the assumed scaling 
the contribution R^ can be bounded by an exponentially vanishing term that is 

( Fi) 

independent on g; all the restrictions on the range of g come from R N , accordingly 
they are the same as in the previous theorems. 

It is clear that with scaling Eq. ( |3.73|) the previous bound Eq. (|3.68|) is too loose 

to prove the convergence: we thus step back to Eq. ( |3.65| ) that is common to both 

theorems and try to get a better bound. By use of asymptotic behavior Eq. (|3.55| ) 

of Disc S(g), of bounds Eq. ( |3.66| ), Eq. ( |3.67| ), extending the integral range to 3' G 

[— oo, 0] and changing the variable as 8' — *> —ft', we find: 

\3\ N+1 r°° 
\R [ x>\ < const. - p~ c / dS'exp 
1 - 8 °" Jo 



-(N+l + c)$ \P', 



(N + 1 + c)p 1 l b 



(3.75) 



$(/?', u) = log/3' + u K ^ . (3.76) 



where 

;i + /3'W b 

Eq. ( p.75| ) can be estimated by the saddle point method. For fixed u > 0, $ behaves 
as u/(3' l / b near 3' = and as u8 ,<ya ~ l ^ b for 8' — ► oo thus at least one saddle point 
(depending on u and others parameters a, b), Pl(u) exists, defined by 

^$(/t«) = 0. (3.77) 

If (depending on a, b) there are more than one such point, the one corresponding to 
smallest $(/^,w) must be chosen. The bound for Rffi is then 



\R N I < const.- — — p c . e 



-AT0(/3i,-u) 



(3.78) 
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(in the last line powers of N have been neglected). The function $*(«) = (tt), u) 
is an increasing function of u, as can be easily seen from (due to Eq. ( |3.77| )) 

^ = d u ®((3:,u) + d^(p:,u)^ = d u $(fc,u) > (3.79) 

(last inequality comes from a direct inspection of Eq. (|3.76| )). Thus if we restrict 
ourselves to u > w*, with defined by 

$*«) = (u*),0 = 0, (3.80) 

the bound on vanishes exponentially. Q.E.D. 

The value of in Theorem 3 should be obtained (for fixed a, b) by numerical 
solution of coupled equations Q3.77j ) and (|3.80|) . Nevertheless it is easy to get a (more 
conservative) a, b independent range of convergence in terms of u. Indeed, if we 
bound $ in Eq. ( |3.76|) with 

ii 

d>(/3» > = log/5' + — - b (3.81) 

and proceed as in Theorem 3, the equations analogous to Eq. (|3.77| ) and Eq. ( |3.80| ) 
can be solved exactly, obtaining the exponential convergence ($(/3^,,m) > 0) for 

u > be' 1 > e _1 . (3.82) 



4 Convergence of the Optimized (Variational) Expansion 

Delta expansion - order dependent mappings diverges if the parameters p or x 
depends on iV improperly, i.e., if they lie outside the convergence range, Eqs. ( |3.28| ), 
( |3.53|) , ( |3.32|) , (|3.54| ) . In the case where x grows with iV faster than N 1 ^ 2 , the behavior 
of the remainder Rn can be estimated by observing that its main contribution comes 
from the integration region near /3' = 1 on T' B . One finds 

S» -(£)'. (4-1) 

(note that x/N 1 ^ 2 — > oo; (3 N — > 1). For the particular case a = 1/2 it reduces to the 
result found earlier IHIl. 
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On the other hand, if x grows too slowly (x/N 



l/2a 



0), the main contribution 



to Rn comes from the region /3' ~ . This contribution can be easily estimated by 
the saddle point approximation and gives: 



(this behavior is an offspring of the standard divergence of perturbative coefficients). 

It is clear that there are no solutions of optimization condition (e.g., dS^/dQ = 0, 
or Sjv+i — 5jv = 0) in these ranges of x = VL/uj. Optimization forces the parameter 
x (or p) to lie within the convergence range. A study of the roots of optimization 
equation (see Appendix C) with regard to the g dependence of p or x, gives a further 
consistency check of this point. 

In other words, we have proven that optimized DE-ODM for anharmonic oscilla- 
tor/double well system converges as long as the coupling constant lies in the appro- 
priate domain. This is important since in many (especially in higher dimensional) 
systems it is a hard task to do higher order calculations. Optimization is crucial to 
get a good approximation from lowest order calculations. 

Of course, optimization condition alone does not usually yield a unique solution 
for x (or p), or a given optimization criterion may not have any real solution at all at 
a given order N. Furthermore, as the number of the roots of optimization equation 
in general grows with N, it is not very practical to focus our attention to a particular 
root and to try to follow it: it may disappear as iV is varied and a new pair of roots 
may appear, etc. JJ. 

In the case of the energy eigenvalues of anharmonic oscillator or double well, a 
general criterion that follows from our analysis is that the convergence rate is fastest 
for the smallest possible Q, in accordance with the earlier suggestions ||. 

In any event, the importance of our results in the precedent sections lies in the fact 
that they do prove the convergence of optimized or variational expansion, without 
however letting ourselves involved into the analysis of complicated behavior of the 
solutions of a particular optimization scheme. We find it remarkable that a variational 
method is proven to converge, even if no variational principle is involved here. 
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5 q 2M oscillators with M > 3 

The convergence proofs presented in Sec. [5] do not apply directly to the case of 
higher anharmonic oscillators (M > 3 in Eq. ( |1.1|) ). For simplicity consider the case 
of real g. All the pre-requisites of Theorem 1 are satisfied: E(g) is known to be 
analytic in the cut g plane |5|, and the conditions 2) and 3) are satisfied with 
p = 1/(M + 1) and b = M — 1. The convergence cannot however be proven in these 
cases because there are no values of the scaling index 7' satisfying b < 7' < a, since 
b = M — 1 > 2 on the one hand, and 1 < a < 2 on the other. 

This does not yet mean that DE-ODM diverges in these cases, since Theorem 1 
gives only sufficient conditions for convergence. Let us re-examine the last condition 
(a < 2) imposed on the conformal transformation. It is needed so that the complex 
part (r#) of the image of the negative real g axis lies outside of the unit circle, i.e., 
\(3'\ > 1, V/3' g r B , p + I. 

Now one knows that E(g) is analytic on the cut g plane, hence E(g(f3)) analytic 
within V surrounded by Ta, T b . If the function \l/(/3,p) = (1 — j3) a E(g) would be 
analytic in a slightly (but sufficiently) larger region of /3, one could consider a con- 
formal transformation with a > 2 (actually one needs a > M — 1), deform (enlarge) 
the integration contour Ta, r# to T' A , T' B until \j3'\ > 1, V/5' G T' B . All the steps of 
the convergence proof would then be valid. Note however that such a deformation 
in the (3 plane involves going onto the higher Riemann sheets in the original g plane. 
It is here that Bender- Wu type singularities potentially make obstructions. Unfortu- 
nately, having no detailed knowledge on the positions of Bender- Wu type singularities 
in the general q 2M oscillators, we are at present unable to answer whether or not DE- 
ODM for q 2M (M > 3) oscillators converges with an appropriate scaling of p (or by 
optimization). 

It is interesting that Pade's method does not converge [0] for q 2M oscillators, for 
M > 3, either. 

6 Green's functions in the one dimensional 4 theory 

As an application of Theorem 1 of Sec. \i.'d[ consider the Euclidean Green's func- 
tions of the anharmonic oscillator system: 

(q(ri)q{r 2 ) ■ ■ ■ q(r n )) 
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- J J Dgg(r 1 )g(r 2 ) • ■•q{r n ) exp 




-,G n ( -^-]UJTi,UJT2, ■ ■ -UJT n ) . (6.1) 



- cW 2 n W 

A simple-minded application of the delta-expansion based on the shift of the quadratic 



term of the Hamiltonian Eq. ( |2.5|) , seems to be problematic, due to the compli- 
cated dependence on cjTj. One may, however, well apply Theorem 1 directly to 
G n (g, {cuTj})!^/^ treating {uJTi} as fixed parameters. 

We first assume the analyticity of G n (g, {uJTi}) in the whole cut plane of g. This 
is quite reasonable, since any given Green function can be expressed as an (infinite) 
sum of terms analytic in the cut plane, by expanding in the energy eigenstate basis. 
Such an assumption was also made by Strharsky [|T5|] to compute successfully the large 
order behavior of the Euclidean Green's functions by use of the dispersion relation. 
The other conditions 2) and 3) of Theorem 1 are satisfied, with p = — n/4; a = 4/3; 
b = 1; c = n/2. 

DE-ODM for G n constructed from its standard perturbation series by the confor- 
mal transformation Eq. Q3.56 ), followed by the Taylor expansion in /3, then converges 
to the correct answer, if p is scaled with the order N as in Eqs. ( 3.28| ), ( p. 29 ) or 
Eqs. (ET53) , (|3733| ). 

The convergence proof above holds pointwise for each fixed value of {o;Tj}. In fact 
a uniform convergence over arbitrary value of {uTi} is desired, because the functional 
form of G n , {uJTi} as the argument, is what one wishes to reproduce. The following 
argument can be given on this point using the result of [0. 

If we apply the same procedure as in ]I| also for G n , we have an expression for the 
remainder similar to Eq. (2.53) of [|TJ. A part (corresponding to above) of the 
remainder is bounded (essentially) by Ci/3 N , where c\ is the first order perturbative 
coefficients for G n . Another part (corresponding to above) is estimated by 

the Euclidean bounce solution for g < 0. Now c\ is computed by using the free 
propagator, (<?(T)g(r / )) = e~ w ! r ~ T ''/(2c<j), thus obviously is bounded for the whole 
range of {uJTi}. On the other hand, {uJTi} dependence of R^ is proportional to (for 
N large) /f^, dr fli v^2sechu;(rj — r ), where r is the position of the bounce center, 
and this is also bounded for arbitrary value of {o;Tj}. This argues strongly for the 
uniform convergence with respect to {u;Tj}. 
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7 Summary and discussion 



Delta expansion - order dependent mappings, is a powerful resummation technique 
which converts a given divergent series, under certain conditions, into a convergent 
sequence of approximants. For its generality, the extreme simplicity of use (in par- 
ticular, combined with an optimization procedure), and for its impressive success in 
simple quantum mechanical systems, it seems to be worthwhile to study carefully the 
mechanism underlying such a method. 

In this paper, the convergence proof for the DE-ODM given previously by the 
present authors for the anharmonic oscillator has been refined and its domain of ap- 
plicability considerably extended. The key step for the improvement was the use of 
Cauchy's theorem, written directly in the complex plane of the conformally trans- 
formed variable, instead of making a conformal transformation in the standard dis- 
persion relation. Use of the knowledge on the analyticity property of E(g) on the full 
Riemann surface was also crucial. 

In particular, this allowed us to show the uniformity of convergence in a wide region 
of the Riemann surface (Fig. 3), including the strong coupling limit of anharmonic 
oscillator (g = oo), and the large g regime of the quantum mechanical double well. 

It is interesting to note that in this last case the standard perturbation series is non 
Borel summable. Together with the example of the zero dimensional model discussed 
Appendix H, our results convincingly demonstrate the fact that the non Borel 



in 



summability of the standard perturbation series is no obstruction to the existence of 
a simple resummation procedure which converges to the exact answer. The possibility 
of reconstruction in non Borel summable cases here hinges upon the known existence 
of a wide region of analyticity of E(g), which is larger than normally required for the 
Borel summability Jl6| . 



As for the variational aspect of the delta expansion - order dependent mappings, 
it is of considerable interest that a variational method can be proven to converge, 
without use of the variational principle. 

A lesson to be learned from the results found here is the fact that, if only con- 
vergence for real finite g is required, quite a wide arbitrariness exits in the choice 
of the conformal transformation, to be used in the construction of DE-ODM. Our 
convergence proof can probably be generalized to conformal transformations of more 
general forms. 
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At the same time, we have clearly marked the limit of validity of the DE-ODM 
based on a simple conformal transformation of g. In particular, the perturbation 
series of the quantum mechanical double well in the regime of weak coupling (where 
tunnelling effects — or instantons — are clearly essential), cannot be resummed by the 
present method. 

Also, in four dimensional field theory models such as Quantum Chromodynam- 
ics, where the analytic structure in the renormalized coupling constant is far more 
complicated ||17|| , the validity of an analogous approach is not obvious. 

Nevertheless, we believe that our general theorems apply to many quantum me- 
chanical and field theory models in dimensions less than four. 

Note Added 

When this paper was being typed the authors were informed that an analogous 
technique was used by H. Kleinert |1| to show the convergence of DE-ODM in the 
strong coupling limit. 
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Appendix A Some properties of the conformal transforma- 
tion 

We study here how the negative g axis is mapped by the inverse of the conformal 
transformation 

9 = P , 13 ^ (A.l) 

for different values of a > 1. The problem is equivalent to the study of solutions of 
the equation 

H(p,s) = (1 - (3) a + s(3 = (A.2) 
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with s = —p/g > 0. 

For s > s c = a a /(a — Eq. ( |A.2| ) has two real solutions that coalesce (one 

coming from (3 = — oo, the other from (3 = 0") for s = s c at f3 — (3 C = — l/(ot — 1). 
These values are found by imposing H((3 C , s c ) = dpH([3 c , s c ) = 0. For < s < s c the 
equation (|A.2| ) has two complex conjugate solutions that start from (3 = j3 c and join 
at (3=1. 

Other complex conjugate solutions converging to (3 = 1 at s = are in general 
also present ,Q but the complex contour that wraps around the cut g plane is mapped 
to the contour T formed by the two above mentioned complex conjugate paths (called 
Tb in Sec. H) pl us the path [f3 c , 0] (r^) of the smaller negative root. 

First of all, solving Eq. ( |A.2| ) in the limit s — > (i.e., (3 ~ 1) one easily finds that 
r# leaves the point (3=1 with angle <p a = ±(1 — l/ot)Tt. 

To study better the shape of T B , set (3 = \(3\e lS and consider \(3\ at a fixed phase S, 
< 5 < 7i, as a function of a. One finds from Eq. ( |A.2| ), after eliminating s, 



mi ■ 7T - 5 / . (a — l)(7r — 5) .. . 

L9 =sm /sin- ^ ^. (A.3) 

a / a 



It is then easy to see that 



\(3\ > 1, if 1 < a < 2; 

< 1, if a > 2, (A.4) 



for V<5, < 5 < TV. 



Appendix B Zero dimensional model 

In this appendix, we consider the delta expansion or, equivalently, the order de- 
pendent mapping with conformal transformation, 

g = p(3j{l-(3)\ (B.l) 

of a simple integral, 



Z(g) = ^= r dqe-o 2 -^ (B.2) 
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It can be shown that there exist 2(fc + 1) roots on the first Riemann sheet of the mapping 



Eq. (A.l), — 7r < Arg (1 — (3) < 7r, in which k is a nonnegative integer such that 2k + 1 < a < 2k + 3. 



For the case of our main interest, 1 < a < 2, therefore only two roots exist on this sheet. 
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(the partition function of 0q theory). Due to the absence of Bender- Wu type singu- 
larities, DE-ODM in this case converges for each g in the whole Riemann surface, in 
the limit, N -> oo, p = C'/N^' -> with 

1 < i < 2, or 

y=l; C" < C^ it = 1.116711531873221 (B.3) 

(In terms of x, x ~ (g/p) 1 ^ 4 ~~ ► °°-)Q 

The proof for (7 lying on higher Riemann sheets in this case requires a non-trivial 
extension of the treatment given in Sec. |[ 

We start with the observation that the analytic continuation of Z(g) is given in 
terms of the parabolic cylinder function by 

z(g) = (^T74 el/(89) ^-v 2 (V V^) (B.4) 

with the strong coupling expansion 

z{g) = ^- 1/4 E(-ir r(n/2 | 1/4) (^ 1/2 )" (B-5) 

^V 71 " n=0 n ' 

(the latter can easily be derived directly from Eq. (p.2j )). 

Following the general procedure described in Sec. [T2] we introduce the function 

9(J3,p)by 

Z(g) = (l-f3) 1 ^,p), (B.6) 

where the Nth. order approximant for Z(g) is to be constructed as Z N (g) = (1 — /5) 1 / 2 • 
{ty(/3, p)} N ; {f(/3)}k is the kth order truncation of the Taylor expansion of /(/?) in 

In order to write Cauchy's theorem in the (3 variable, the analytic property of ^ 
must be studied first. A good knowledge of it can be deduced from the strong coupling 
expansion Eq. ( |B.5|) . Indeed, it is easily shown that this expansion converges for any 
g 7^ 0, implying that the function \1/ defined above is analytic everywhere, except 
at the quartic branch points at — and at (3 — 00. (j3 = is also an essential 
singularity, see Eq. flB.ll|) below.) 

12 Previously convergence was proved only pointwise for real g || [lj and for a particular phase 
of g, Argg = 2ir (corresponding to the double well) |2|. 
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The image of the cut g plane is the interior of a unit circle (minus the negative 
real axis) in the first (3 Riemann sheet and it is tempting to take its boundary as 
the integration contour r = Ta + ^ b and make a "small" enlargement of the contour 
around (3=1, keeping the rest of the contour unchanged, as was done in Sec. EJ72 



It turns out that this procedure does not work here. (The convergence factor 1/(1 + 
rocos0 a + 0(tq)) in Eq. (|3.72 ) is lost in the case a = 2.) On the other hand, the 
analyticity domain of \I/ is much wider here, so that the contour can be enlarged more 
substantially. Such enlargement is indeed possible but involves a delicate issue; it can 
be dealt with as follows. 

We write the remainder Rn(P) = Z(g) — Z N (g) as: 

BwOT = (1 -,)V V+1 ^ +r J_?M_, (B . 7) 

where T' B will be taken as a circle with the radius R > 1 and T' A as a contour which 
wraps around a part of the negative real axis —R < /3' < (see Fig. 4). 

The radius R will be taken dependent on the order N as 

R=l + dN~ K , (d>0; i/2 > « > 0). (B.8) 

In deriving Eq. (|B.7| ), it was implicitly assumed that \(3\ < R. This indeed holds 
for large iV with this choice of R. ((3 ~ 1 — (p/g) 1 ^ 2 as p — ► 0.) (For the subtleties 
involved in writing Cauchy's theorem on (3 Riemann surface and for its elucidation, 
see the comment before Eq. Q3.40 ) in the main text.) 

On T^, \f3'\ = R and — /3| > R — hence the contribution to Rn(P) is 
bounded by 

Kw ,i s|1 - fl v.(^«5gi 

When (3' is on the contour T' B , the corresponding g' = p/3'/(l — f3') 2 circles around 

the origin g' = on the second Riemann sheet, with the radius, 

, ii R R / \ 

M = '^^, IS 3 TF (B.10) 

With p scaling as p = C'/N' 7 , and R as Eq. ( |B.8| ), \g'\ is small for large. We can 
therefore use the asymptotic expansion near the origin,^ 

^ J_ g r(2n + l/2) r ± .^ VW ] „ (Rn) 
n! L J 



13 The plus sign is for —5ir/2 < Argg < —tt/2 and the minus sign is for ir/2 < Argg < 5tt/2; the 
terms with the exponential factor is absent in other regions of g [^(| Eq. (9.246)]. 
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to get the following bound on 



W0I<777^Lt^1> 



(R-iy/ 2 



(R 



ApR 



[l + 0(pR/(R-lf)). 



[B.12) 



We get 



< 



V2\l-P\ 1/2 



(\(3\\ N+1 


[(^-1) 2 1 


{-r) exp 


ApR 



~ exp I -dN 1 -" + 



AC 1 ' 



■N 1 



'-2k 



(B.13) 



Note that for < 7' < 2, it is always possible to choose n in such a way that 

0<Y-2k<1-k. (B.14) 



With the first factor dominating, one gets \R 



N 



as iV — > 00. 



On T^, — i? < /?' < 0. The value of g' corresponding to it moves along — p/A < 
g' < 0. Since p <C 1 for iV large, the asymptotic form 



|Im*(/3')| 



can be used in 



|ImZ(^')| 1 e 1 /^') 
|1 -p'] 1 / 2 ~ 72|1-/5'| 1/2 



1 + W)), 



1/2/qiV+l 



Therefore, using \(3' — (3\ > \(3\ and |1 — f3'\ > 1, we get 



(B.15) 



(B.16) 



< 



V2 l 



p\ 1/2 m N 



dp' 



7T 



^+W)(i + (l/iV)), 



(B.17) 



where $(/?') = In/?' + u(l + /3') 2 /(3', u = 1/(4(JV + l)p) and in the last step we have 
changed the variable /?' — > — /3' and extended the upper limit of integration to 00. 
The right hand side of Eq. (|B.17|) can be evaluated by the saddle point method for iV 
large, as was done in Sec. |^, Sec. Combining the bound for R^ thus obtained 
with that for R V N , we get the announced result. 

14 We have taken the dominant term for Re g' > 0. When Re g' < 0, i s bounded by a power 

of \g'\ and the contribution is harmless for the convergence. The existence of such a large factor 
near the origin on the second Riemann sheet is understood as g = e 2nl \g\ corresponds to double well 
model. 
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An interesting application of our result is the case of the "double well" model in 
zero dimension: 

z{g) = -j= j™Ji^- 9qA = T^xr^^d- 1 /^), (B.18) 



with g real positive. From the known properties of the parabolic cylinder function one 
finds that Z(g) = {e 2m ) l ^Z{e 2m g) (an analogue of the Symanzik's scaling relation). 
Consequently, the (scaled) delta expansion of the double well model (whose standard 
perturbation series is non Borel summable) converges to the exact answer for any real 



value of g 7^ 0, in accordance with an earlier result [19 



As pointed out in 0, the case Argg = 7r can also be related to Z(g). The crucial 
relation is 

Z{g) = V2e 1/{4g) ReZ{e m g). (B.19) 

Again according to our result, the exact answer for Z(g) can be correctly reconstructed 
systematically by using DE-ODM for Z(e m g), whose standard perturbation series 
Z(e ni g) ~ (1/v 7 ^) En=o T ( 2n + l/2)g n /n\ is non Borel summable. 



Appendix C Optimization condition on p 

In this appendix, we show that either of optimization conditions, the principle of 
minimal sensitivity (PMS) or the fastest apparent convergence (FAC), leads a simple 
equation for p which is independent on g. This fact, for the particular case with 
a = l/2 and a = (M + l)/2, has recently been noted by [pi}] . 

By PMS, one imposes dS^ DM ^ / dx = (or equivalently dS^ DM ^ / dp = 0) order by 
order. From the definition of DE-ODM, we know 

N 



1 

\ L P) n =0 

N 

£ c n g n (l - P) an -° {(1 - P)- an+a } N _ n , (C.l) 



n=0 



where {/(/3)}fe is a kth order truncation of the Taylor series on (3. With a use of a 
relation 

I [d - m' {a - = - £ ^p(i - /rv, (c.2) 



34 



we have (note the variation with respect to x is taken while g is kept fixed) 

^c(odm) W (3 N ^ T(N+(a-l)n-a + l) n 

dx N dx(l-PY^^ Cn {N-n)\Y{an-a) P ' 1 j 

PMS therefore leads an equation 

^ r 'l +(a M"°l" ^a (C4) 

(N -ny.r(an-a) 

On the other hand, it is easy to see that FAC, S'^ )DM ' ) — S'^ D 1 M ' ) = 0, leads a similar 
equation 

£o [N -ny.r(an-a) 

Optimization condition thus leads to an order by order determination (albeit 
among all possible, at most iVple, roots) of p which depends only on N but inde- 
pendent of g. On the other hand, our convergence proof is valid if the parameter p 
is appropriately scaled with N, but independently of g. See Eq. ( |3.52| ) or Eq. ( |3.53| ). 



These facts are consistent with each other: as noticed in Sec. f|, optimization forces 
x or p to lie within the range Eq. ( |3.52| ) or Eq. Q3.53 ) for which DE-ODM converges, 



hence selects out p which grows appropriately with iV but in a g independent manner. 



Appendix D Critical exponents in 0| theory 

In this appendix, we apply DE-ODM to the renormalization group functions of 
03 theory and compute various critical exponents. This problem was already studied 
in [0], but our approach here differs in the following aspect: In 0, a renormalization 
group argument was used to relate first the analytic structure in g Q around go = oo (go 
is the unrenormalized coupling constant) to the behavior of the beta function W(g) 
near the infrared fixed point g ~ g*. The parameter of the conformal transformation 
is accordingly taken to be a = 1/oj (u being the slope of the beta function at the 
fixed point), in order to render W (g(go(f3))) analytic around (3 = 1. The input value 
of uj must somehow be guessed by e.g., other approximate methods for computing 
W(g). Although quite encouraging results have been found 0, the fact that one must 
assume an input value of uj (which is one of the results being sought for) to start the 
calculation with, is a somewhat disturbing feature of the approach . 
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On the other hand, as is seen in Sec. |3.1| , the requirement of analyticity near 
(3 = 1 is not really necessary, if one is only interested in the convergence for finite 
real positive g. Also, quite a wide range of conformal transformations can be used in 
such a case. In this spirit one may a priori try a simpler conformal transformation, 
directly to the renormalized coupling constant g. 

Here we choose a = 3/2 and a = 1/2, the same value as in the case of the 
anharmonic oscillator. Assuming the renormalization group functions to have the 
analyticity on the cut g plane (which is not proven to our knowledge), one concludes 
that DE-ODM gives a convergent sequence for true functions.! 15 ] 

We thus proceed as follows: take a perturbative series of, say, the beta function 

oo 

W(g)~Y. c n9 n , (D.l) 

n=0 



for which the first six non-trivial coefficients have been calculated | 25(| . Then ODM 
with a = 3/2 and o = 1/2 gives a new series 

„, i < £ » 1 ^" r(3n/2 + * - 1/2) / 1\* 

WMs,x) = Sc., g t , r(3n/2 _ 1/2) (1 - ~ 2 ) , (D.2) 



where p and /? have been parameterized in terms of x as in Eq. ( |3.2| ) . At lower orders 
we expect that an optimized determination of xn gives a more accurate result than 
with a particular scaling of xn with N. We thus employ here the principle of minimal 
sensitivity, i.e., dW^/dx^ = 0, for each value of g, by choosing the root closest to 
x N = I. 

In Fig. 5, we plot the original beta function Eq. ( p.l|) up to 0(g 7 ) p5 and the 



function Wi{g , x<r{g)) obtained in this way. A clear improvement and the appearance 
of an infrared fixed point can be seen. 

We applied the same method for perturbative series of various renormalization 
group functions |2"5|] . The resulting values of the critical exponents, as a function of 
the order N, are summarized in Table The first two columns of Table, g* and lu, 
are the result of the method explained above. In the last two columns, we have used 

15 It is known that the perturbation series for the renormalization group functions is asymptotic 
for Re g > p3| , and the known large order behavior of those functions |M is consistent with the 



assumed form of the imaginary part Eq. (3.55) with 6 = 1. 
16 For the definition of various critical exponents, see [Em. 
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g* = 1.43645, the result of the seventh order resummation. The blank means there is 
no solution of the optimization condition near x^ ~ 1. 

These results can be compared with those obtained by the Borel summation 



method |pq| : 

g* = 1.416 ±0.005, cu = 0.79 ± 0.03, 

rj = 0.031 ±0.004, 7 = 1.241 ± 0.0020. (D.3) 

Considering the simplicity of our approach (in particular, the uniform convergence 
for g G [0, 00] is not guaranteed), the agreement is rather surprising and suggests the 
correctness of our assumptions. 



Appendix E Strong coupling expansion coefficients from stan- 
dard perturbation series 



As shown in Sec. |3.2| DE-ODM for the anharmonic oscillator converges even in 
the infinite coupling constant limit. This fact can be used to reconstruct the strong 
coupling expansion Eq. ( |1.4j ) from the perturbation series coefficients c n [12]. Here 



we derive such a formula, which is more general than that given in |L2] and whose 
derivation is also somewhat simpler. 

It is interesting that a quite similar formula had been proposed earlier [IT], fol- 
lowing a completely different approach. 

The iVth order DE-ODM approximant is given by |l j 



N 



N-n 



T(3n/2 + A;- 1/2) 
n=o - k =o fc!r(3n/2-l/2) 

Substituting x = g 1/3 y N (see Eqs. (|333|) , (gjp ) into Eq. 



(E.l) 



S 



n = 9 



1/3 



N N-n k 

£££> 

n=0 k=0 1=0 



T(3n/2 + k — 1/2) 



A\), we find 
(-1)' 



T(3n/2-l/2)(A;-/)!^+ 2i - 1 l\ 



-{9 



2/3\l 



(E.2) 



A comparison with Eq. ( |1.4| ) then gives for the /th strong coupling expansion coeffi- 
cients di, 



d. 



lim 



I) 1 ^ N d^ n r(3n/2 + k + l — 1/2) 1 



n=0 k=0 



k\T(3n/2 - 1/2) 



Dn 



(E.3) 
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(-l)'^f r(iV + n/2 + l/2) 



^-°o /! n (A^-/-n)!r(3n/2- l/2)(3n/2 + /- l/2)^+2«-i' 

where yjy must be scaled appropriately with N (y^ = ciV 7 ; 1/3 < 7 < 1/2 with 
arbitrary positive c, or y^ = ciV 1//3 ; c > c* = 0.570875 •• •), or simply determined 
order by order by optimization of Sn or d\. 

Delta expansion for the double well potential is obtained simply by changing 
(3 — > 1 + 1/a; 2 From the above derivation, it is obvious that e^ DW ^ = (— l) l d\ AUO \ 
which is of course a correct result, in view of Eq. ( jl.3| ). 

Finally, as a further check, the known strong coupling expansion ( |B.5|) in the zero 
dimensional model discussed in [Appendix B| can exactly be reproduced following the 
same procedure. 
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Table Caption 

1. Determination of various critical exponents in 0g theory by DE-ODM. 
Figure Captions 

1. The domain of analyticity T> in {3 plane, corresponding to the cut g plane. Its 

boundary is the curve Ta which wraps around the negative real (3 axis, and Tb, 
the upper and lower complex sections. 

2. Enlargement of the contour Tb to T' B , so as to include part of the circle C. The 

crosses represent the (schematic) position of Bender- Wu singularities. 

3. The inverse image of the integration contour T' B (solid line) and the boundary of 

the region Eq. (|3.49| ) (dashed line) on the g Riemann surface. Straight sections 
of the contour image are actually on the axis, Argg = ±7r, ±2tt, ±4tt. 

4. The contour used in the convergence proof of the zero dimensional model. 

5. Standard perturbation result for the beta function W(g) (dotted line) and DE- 

ODM improved one (solid line). 
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N g* lo r] 7 

2 1.22847 

3 1.77550 0.548573 0.025686 

4 1.24137 

5 1.48560 0.699459 0.030055 

6 1.24366 

7 1.43645 0.751379 

Table 1: Determination of various critical exponents in 0g theory by DE-ODM. 
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